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1 Introduction 



The goal of this paper is to present a quantum adiabatic theorem that is general enough to 
apply to situations in which Hamiltonians have resonances instead of bound states or if just 
an approximation to a bound state is known. The hypotheses of our main result, Theorem 
11.11 do not specifically mention resonances, so we demonstrate how one applies the result 
by considering the specific situation of shape resonances. We plan to apply our theorem to 
other resonance situations in the future. 

Our application to shape resonances has considerable overlap with the work of Abou- 
Salem and Frohlich pj, although many of the details are quite different. In some instances, 
we obtain sharper estimates. 

The adiabatic theorem of quantum mechanics describes the long time behavior of so- 
lutions to the time-dependent Schrodinger equation when the Hamiltonian generating the 
evolution depends slowly on time. The theorem relates these solutions to spectral informa- 
tion of the instantaneous Hamiltonian. 

The traditional quantum adiabatic theorem applies to Hamiltonians that have an eigen- 
value which is separated from the rest of the spectrum by a gap. Some more recent versions 
do not require the gap condition. All one really needs for the adiabatic theorem is a spec- 
tral projection for the Hamiltonian that depends smoothly on time. This allows situations 
where an eigenvalue is embedded in the absolutely continuous spectrum of the Hamiltonian. 
Since embedded eigenvalues are intrinsically unstable, they usually become resonances once 
the system is perturbed. It is intuitively clear that on the time scales which are shorter 
than the resonance lifetime, there should not be much of the difference between a proper 
bound state and the corresponding resonance. This intuition leads to the question whether 
an adiabatic theorem holds if there is a nearly spectral projection for the Hamiltonian that 
depends smoothly on time. The main abstract theorem of this paper. Theorem II. H provides 
an affirmative answer to this question. 

The paper is organized as follows: We state our abstract result in Section 1.1. We 
then describe its application to shape resonances in Section 1.2. Section 2 contains the 
proofs. Technical details are collected in the Appendix. 
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1.1 The Extended Adiabatic Theorem 

To state an adiabatic theorem precisely, it is convenient to replace the physical time t by the 
rescaled time s = et. One is then concerned with the solution of the initial value problem 

teij,{s) = H{s)^,{s), with Ve(0) = i^o, (1.1) 

for small values of e. The Hamiltonian H{s) is self-adjoint for each s and depends sufficiently 
smoothly on s in an appropriate sense, and ip^ takes values in the Hilbert space. We shall 
be more specific about what we mean by smoothness below. Typically, s is kept in a fixed 
interval, so that the physical time t belongs to an interval of length 0(e^^). 



Our main result. Theorem 11.11 hinges on three assumptions given below. The first is a 
static condition. The second is a dynamic condition imposed on the family H{s). The third 
controls the relative boundness of the rate at which H{s) changes. 

Definition We say that a projection P{s) is nearly spectral for H{s) if it is self-adjoint 
and 

II {H{s) - E{s)) P{s) II < 5/2, (1.2) 
for all s G [0, 1], where 5 is a small parameter. 

Definition We say that a projection P{s) is smoothly nearly spectral for H{s) if it is 
self-adjoint, P{0) is nearly spectral for H{0), and 



^^[{H{s) - E{s)) P{s] 
for all s G [0, 1]. 



< 5/2, (1.3) 



Remarks 

1. If P{s) is smoothly nearly spectral, then it is nearly spectral. 

2. If P{s) is the spectral projection for energy E{s), then it is smoothly nearly spectral 
with 6 = 0. 

Assumption 1. We assume that there exists either a smoothly nearly spectral or a nearly 
spectral projection P{s) for H{s). 
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Remarks 



1. As we have already mentioned, it is reasonable to expect adiabatic behavior of the 
system whenever e is small, but 5. 

2. From the Weyl Criterion (Theorem VII. 12 of [I2]), a non-trivial nearly spectral pro- 
jection exists for any point E in the spectrum of a self-adjoint operator. However, we 
shall impose further dynamical assumptions on P{s) below. In general, the dynamical 
assumptions limit the set of suitable E{s) to eigenvalues or resonances. 

3. The notion of a nearly spectral projection is also related to the ideas of Spectral 
Concentration. See, e.g., Section XII. 5 of [H]. 

Assumption [1] is a static condition imposed on the family H{s). Our results require a 
dynamic hypothesis as well: We let Qa be a smoothed characteristic function that takes the 
value 1 at 0. More precisely, we assume 

ga{x) = 1 if |x| < a/2, ga{x) = if |x| > a, and ga G C^i^IR). (1.4) 
Assumption 2. There exists a G (0, 1) such that 



uniformly for s E [0, 1]. 
Remarks 

1. Here S' is another small parameter, while a should be thought of as an auxiliary 
tuning mechanism, which eventually could be optimized. In our application to shape 
resonances 6' is roughly of the same order of magnitude as 6. 

2. Intuitively, condition (11. 5p quantifies the rate at which the wave function "leaks" from 
the range of P{s) to energetically close states. One can also think of Assumption [2] 
as requiring a bound on the spectral density in the vicinity of the resonance or bound 
state. Indeed, if there are no bound states with energies close to that of the resonance, 
the expression in (II. 5p tends to zero as a tends to zero. The second possibility is that 
there are bound states nearby {e.g., pure point spectrum), but their overlap with P{s) 
is small. The latter occurs in the Anderson localization problem. 
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3. We remarked earlier that non-trivial nearly spectral projections exist for any self- 
adjoint operator. However, if the operator depends on a parameter s, it may be 
impossible to construct a corresponding family of nearly spectral projection P{s) that 
satisfies Assumption |2l Theorem 11.21 shows that both Assumptions [1] and [2] can be 
satisfied for shape resonances. 



Finally, because we are dealing with unbounded perturbations, we impose the following 
technical requirement: 

Assumption 3. We assume that the operators H{s) are self-adjoint on a common domain 
V, and that the derivatives H{s) and H{s) exist as operators from T> to the Hilbert space "H. 
Furthermore, we assume there exists C , such that for all s G [0, 1], 



H{s) (His) 



< C 



and 



His) (His) - i) 



-1 



< C. 



Here and throughout the text, C stands for a generic constant. 

Our main abstract theorem is the following: 

Theorem 1.1. Suppose P{s) is a nearly spectral projection for the Hamiltonian H{s) that 
satisfies Assumptions\^ -0 Then there exists a unitary propagator for M.l\) . 
Let V^,(0) G Range P(0). 

If P{s) is nearly spectral, then for all s G [0, 1], 



dist{7/',(s), RangeP(s)} < 2 6' + 2 \\P{s)\\ \\\ga\h6 + K^e + 6/e. 



1.6) 



If P{s) is smoothly nearly spectral, then for all s G [0, 1], 



C6 



dist{7/',(s), RangeP(s)} < 2 6' + 2 \\P{s)\\ \\\ga\hS + ^ + C K^e + 6/e. (1.7) 



The a-dependent constants Ka and Ka in these expressions are given in (\2.8\\ (respectively 
f l2.10p ) below. The norm \\\ga\\\3 is one of the norms described in 16] and Appendix B of [11]. 
These norms have the form 

\\\9a\\U2 = W^Wk-n-l^ ^^^^ 11/11' = / a+Xy/'\fix)\dx, 

k=0 



where the C depends only on n. 



Remarks 



1. If -P(s) is a spectral projection, we can take 6 = 0, and ( II. 7p corresponds to the 
(slightly improved) adiabatic theorem of [5]. Even in the case of a bound state it 
is often technically simpler to construct an approximant P{s) for the bound state 
projection rather than the exact eigenprojection. The error in the approximation then 
contains the parameter 6. 

2. If P{s) corresponds to an eigenprojection for an isolated eigenvalue of H{s), then it 
is typically relatively easy to verify its smoothness (with respect to the s variable), 
as it is "inherited" from the smoothness of H{s). Otherwise, checking this is usually 
highly non-trivial (c./. |9] where this task is carried for the ground state of the atom 
in a QED picture). One of the main obstacles in the implementation of the quantum 
adiabatic algorithm [7j is controlling the norms of the derivatives of P{s). From this 
point of view, (11. 7p is a better result than (11.61) . as it only requires a bound on the first 
derivative of P{s). We note that even if P{s) is only nearly spectral, one can get a 
bound in the adiabatic theorem that depends only on the norm of ||P(s)||, but not 
on P{s). That can be achieved using a mollifier argument. In our application to shape 
resonances, we have the good control on the smoothness of P{s). (See Lemma [2. 3[ ) 

1.2 The Application to Shape Resonances 

Theorem 11.11 applies to Schrodinger operators with shape resonances. The analysis is not 
much more complicated for certain magnetic fields, so we include them. We consider 
Schrodinger operators H{s) := 'Pa 'Pa + ^(s) on iR'^, where = — i /i V — A(x). 

Assumption 4. We assume the components of the vector potential A and their first 
derivatives diA for i = l,...,d are bounded on IR'^. 

Shape resonances (see, e.g. [1]) are resonances of H{s) that arise because the particle 
can be confined to a region of space that is bounded by a classically forbidden region. If the 
resonance has energy near E, then we define the classically forbidden region to be 

J := {xelK^ : V{x) > E + b} 
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for some 6 > 0. One usually assumes that J separates IR into a bounded interior and an 
unbounded exterior component. The intuition is that the particle spends a long time in the 
interior component, but can eventually tunnel to the exterior component. 

We examine this situation where the energy E{s), potential V{x^ s), and classically 
forbidden region J{s) all depend on s. 

For simplicity, we assume that for an appropriate value of E{s), J{s) separates iR"' into 
an exterior region 0{s) and a single connected interior region /(s), so that 

= J{s) U 0{s) U /(.). 

(Se( 




Figure 1. A two dimensional potential that has shape resonances. 




Figure 2. Contour plot for the same potential showing the classically allowed regions (red, 
solid curves) and classically forbidden region (blue, dotted curves) for a particular energy. 
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We assume that H{s) has the following properties for each s G [0, 1]: 
Assumption 5. 

1. For every s G [0, 1], we assume H{s) is a self-adjoint operator on an s-independent 
domain T>. 

2. For every s E [0, 1], there exist c > and an open set C JR!^, such that 

dist {0(s), > c, and dist{/(s), Q%s)} > c. 

Moreover, the Friedrichs extenstion Hq{s) of the Dirichlet restriction of H{s) to Q{s) 
is bounded from below. 

3. The operator Hq{s) has a discrete eigenvalue E{s) that depends smoothly on s. 

4- E{s) is separated by a spectral gap A(s) > A > from the rest of the spectrum of 
Hq{s)- For convenience, we assume A < 1. 

5. The potential V{s) satisfies \\ Vj{s) \\ < C , 

V{s) iHnis)-t)-' < C, and || t>(s) (Hnis) - 1)'' \\ < C, 

where Vj stands for the Dirichlet restriction to the set J. 
We further assume that supp V{s) C /(s) . 

Remarks 

1. In typical situations, part 4 of this assumption forces A to be 0{h). 

2. Assumption [5] implies a "Combes-Thomas" estimate for The (improved) 
Combes-Thomas estimate [8] is usually stated for the operators acting on the whole 
space IR'^. The extension to the sub-domain case is presented as Theorem 13.41 in the 
Appendix. 

The following theorem and its corollary are our main results concerning shape resonances. 

Theorem 1.2. Suppose that Assumptionl5\is satisfied. Then there exists a family P{s) of 
projections that satisfies AssumptionsU\-{Mfor energy E{s). The values of the corresponding 
parameters are 

S = Che-'^/\ a = A/2, and 6' = C h e"^/^' IWgaWU/ A , (1.8) 



where rj > 0. 

Combining Theorems 11.11 and 11.21 we obtain 

Corollary 1.3. In the context of a shape resonance, the adiabatic theorem holds for 

e > h-^e-'''^, where > 0. The error term in the adiabatic theorem is bounded by a 

constant times 

Proof The only nontrivial part of the proof of the corollary is estimating factors in (11.61) 
that come from P{s) and -P(s), some of which occur in Ka- However, P{s) from Theorem 
11.21 is constructed in Section [2^ from Pq(s). 

We begin by estimating derivatives of -Pn(s) by writing 

Pn{s) = - -L / {H^{s)-z)-^ dz. 

J\z-E{s)\=Il/2 

The length of the contour is tt A, while self-adjointness and the gap condition assure that 
the norm of the integrand is bounded by 2/ A. 
The derivative is 

Pn{s) = - ^ (f {Hn{s) - z)'' V{s) {Hn{s) - zy' dz. 

J\z-E{s)\=A/2 

By Assumption and the first resolvent formula, the norm of the integrand is bounded by 
C/A\ so \\Pn{s)\\<Cn/A. Similarly ||Pn(s)|| < Ctt/A^. 

The corollary is an immediate consequence of these estimates and Lemma 12. 3[ H 

Remarks 

1. As remarked earlier, A is typically bounded above and below by constants times h. 
When this is the case, the error in the adiabatic theorem is bounded by a constant 
times 

for small H. 

2. There are situations in which eigenvalues of Hq can cross as h is decreased. In those 
situations, we obtain no error estimate. 
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3. As we have already commented, the Weyl Criterion guarantees existence of nearly 
spectral projections at any point of the spectrum of a quantum Hamiltonian. However, 
for general Schrodinger operators, the ranges of those projections contain only functions 
are very delocalized in space. The nearly spectral projections we construct for shape 
resonances are localized in a region that can be chosen independent of 6 and h. 

4. There are numerous definitions of resonances for a quantum mechanical Hamiltonian 
H. One is as follows: Since H is self-adjoint, the complex valued function 

f{z) := {H-z)-'^P) 

is analytic in the upper half-plane for any vector ip. For ip in some appropriate set, / 
has an analytic continuation / to some portion of the lower half-plane. A pole of / at 
E — iT corresponds to a resonance near energy E with lifetime proprotional to l/F. 
This definition is often not practically tractable. 

2 Proofs 

2 . 1 Proof of Theorem [131 

Using Assumption [3l we obtain existence of the unitary propagator U^{s) by applying The- 
orem X.70 of [13]. 

Without loss of generality we can assume that E{s) = throughout the proof. Indeed, 
the dynamics generated by H{s) differs from the one generated by H{s) — E{s) only by 
the (dynamical) phase exp [ — i Jq E{t) dt/e). 

The first step in many proofs of the adiabatic theorem {e.g., [3]) is the construction of 
the so-called adiabatic evolution. This is a unitary family Ua{s), such that 

P{S) = Ua{s) P(0) Ua{sy. (2.1) 

The second step of the proof is to verify that Ua stays close to the true evolution f/^, deter- 
mined by 

i e U,{s) = H{s) U,{s), with U,{0) = I. (2.2) 

We follow this strategy with one modification: Since P{s) is not a spectral projection, it 
is hard to construct an evolution that satisfies (12. ip and is close to Ue. Instead, we construct 
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a nearly adiabatic evolution Un{s) and replace (12 .ip by 

II P{s) Un{s) - Un{s) P(0) II = II P{s) - f/„(s) P(0) V^{s)* II < Sje. (2.3) 

With this modification, the second step is the same as in the traditional adiabatic theorem. 
Specifically, we let Vni^s) be the solution to the initial value problem 

%e\Jn{s) = Hn{s)Un{s), with f/„(0) = /, (2.4) 

where the generator is Hn{s) = H{s) + i e [P{s), P{s)]. Existence of Un{s) is guaranteed 
by Theorem X.70 of [T3] . 

In order to check that the estimate (12. 3p holds, we compute 



^ (^UnisY Pis) Unis)^ = Unis)* [H (s) , P(s)] f/„(s), (2.5) 



where we have used 

P{s) = P{s) P{s) + P{s) P{s) and P{s) P{s) P{s) = 0. 

The condition (11.20 and the unitarity of ?7„(s) guarantee that the right hand side of 
(12. 5 p is bounded in norm by 6/e. We integrate both sides of (12.51) and use the fundamental 
theorem of calculus to see that 

II UnisT Pis) Un{s) - P(0) II < (5/e, (2.6) 

for all s G [0, 1]. The estimate (12. 3p follows from the unitarity of ?7„. 

Next we show that the physical evolution stays close to Un for all s G [0, 1]. 
We claim that for a nearly spectral projection, 

\\Un{s) - U,{S)\\ = \\U,{sy Un{s) - I \\ 

< 2 6' + 2 \\P{s)\\ IWgalhS + Kae, (2.7) 

where 

Ka := 2 \\\ga\\\3 max ||P(r)||2 + 2 |||^„|||3 max ||P(r)|| 

r r 

+ C (\\\ga\h max ||P(r)|| + |||(?,|||4 max ||P(r)||) . (2.8) 

\ r r / 
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For a smoothly nearly spectral projection, we can improve this bound to 

||C/„(s) - U,{s)\\ < 2 5' + 2 \\P{s)\\ WlQalhS + C5/a^ + CK^e, (2.9) 

where 

Ka := max||||<7a|||5 + lll^allk + ^(l +11^(^)11)}- 

To prove these estimates, we note that the unitary operator Q,{s) :— Ue{s)* Un{s) sat- 
isfies the differential equation 

= K(s) n(s), where K(s) := C/,(s)* [P(s), P(s)] Ue(s). 

We cast this in its integral form 

n{s) - I = f K{r) n{r) dr. (2.10) 
Jo 

To show that the right hand side is small, we use the following lemmas that we prove below: 

Lemma 2.1. For a nearly spectral projection P{s), there exist a pair X{s) = -^i(>s) and 
Y{s) = Yi{s), such that 

[P{s), P{s)] = [X,{s), H{s)] + Y,{s), (2.11) 

where 

\\MS)\\ < lll^a|||3 ||P(«)||, (2.12) 

\\Ms)\\ < C (||k|||3||P(^)|| + ||kll|4||P(^)||), (2.13) 

and 

m{s)\\ < 2 6' + 2\\P{s)\\ WMUS. (2.14) 



Lemma 2.2. For a smooth nearly spectral projection P{s), there exist a pair X{s) — X2{s) 
and Y{s) = 1^2 (s), such that 

[P{s), P{s)] = [X2{s), H{s)] + Y2{s), (2.15) 
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where 



\\X,is)\\ < C/a\ 

\\X2m < C (\\\ga\h + lll^?a|||4 \\P{s) 



and 



\\Y2{s)\\ < 2 6' + 2 \\P{s)\\ \\\g,\\\,6 + C5/a' 



(2.16) 
(2.17) 

(2.18) 



If we substitute the representation (12. lip (respectively (I2.15p ) into the right hand side of 
(I2.10p we observe that 



K{r) VL{r) dr 










+ 



U,{ry [X{r), H{r)] U,{r) Vt{r) dr 



(2.19) 



However, 



d 



U^ir)* [X(r), H{r)] f/,(r) = ^ e ^ ( U^ir)* X{r) U,{r) ) - z e U^ir)* X{r) U,{r). 



So, the second contribution in fl2.19p can be integrated by parts: 

f/e(r)* [X(r), H{r)] U,{r) n{r) dr = -it [ U,{ry X{r) U,{r) fi(r) dr 



+ I 6 U,{ry X{r) U,{r) Vt{r) 



- z e / U,{ry X(r) U,{r) l](r) dr. 
'0 



The norm of the first term on the right hand side is bounded by e max^ ||X(s)||. The norm 
of the second term is bounded by 2e max^ ||X(s)||. Finally, since f2(r) is unitary, CL{r) = 
K{r)n{r), and \\K{r)\\ < 2 || P(r) ||, we can bound the last term by 2 e maxJ|X(s) || ||P(s)||. 
Combining these estimates, we get fl2.7p and 02.91) . 

We now let ipo G RanP(O) be a unit vector and set ipe{s) = U^{s) tpQ. Then using 
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O and ([21]), we see that 

Ueis) - Pis)U^is)MO)\\ 

< WMS) - Un{s)tljM\\ + II f/n(s) P(0) ^,(0) - P{s)Un{s)tljM\\ 

< \\Un{s) - U,{s)\\ + 6/e. 

The theorem now follows from (12. 7p ( (12.91) respectively). H 

Proof of Lemmas 12 . II and 12 . 21 The operator valued function ga{H{s)) admits the Helffer- 
Sjostrand representation 

gaiH{s)) = [ {d,ga)iz) Rsiz) dzdz, (2.20) 
where Rs{z) := {H{s) — z) ^ , Qaiz) is supported in the disc \z\ < a, and 

\d,ga{z)\ llmz]-''-' dzdz < \\\ga\\\n+2 ■ (2.21) 



c 



See e.g. [6j or Appendix B of [H]. The norm here is 



lll^allln+2 = H'\_^_,, With 11/11, = / \f{x)\ dx, 

I. n J —OO 



k=0 

where the C depends only on n. 

Our first goal is to show that (II. 2p implies 



[ga{H{s)) - /) P(S)|| < |||(7a|||3 <^/2. (2.22) 

To prove this, we recall that we have taken E{s) = and note that for any z E C, inequality 
( II. 2p implies 



II (His) - z) Pis) + z Pis) II < 5/2. 

Thus, 

- Pis) + iHis) - z)-^ Pis) < 



2 I Imz P 

Substituting this relation into (I2.20p and using the bound (I2.2ip for the error term, we get 
II gaiO) Pis) - gaiHis)) Pis) \\ < \\\ga\h 5/2. 
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Since ga{0) = 1, inequality fl2.22p follows. 

We now define Xi{s) in (l2llD and X2{s) in (12J[5|) to be 

= - / {d,ga)iz) R,iz) P{s) R,{z) dzdz, (2.23) 

and 

X^is) = {l-ga{H{s))) Rl{0)H{s)P{s) + h.c. (2.24) 

Since 

^Rs{z) = - Rs{z) H{s) Rs{z), 
as 

it follows from f l2.2ip and Assumption |3] that Xi{s) satisfies the bounds fl2.12p and fl2.13p . 
Since 

\\{l-g,{H{s))) RliO) {H{s)+z)\\ < C/a\ (2.25) 



and 



Rs{-i) H{s 



< C 



by Assumption [21 ^2(5) satisfies the bound fl2.16p . 

To get the bound f l2.17p we first note that by partial fractions, 

\ Rs{z) = Rs{Of Rsiz) + - Rs{Of + \ Rs{0). 
z^ z z^ 

Second, since [a — z)~^ is the resolvent of multiplication by the constant a, 

9'M = [ ^4^dzdz = [ ^^^dzd-z 

da Jc {z-a) [z - ay 

When a = 0, this is zero. Using these facts and (7a(0) = 1, we can write 

{l-ga{H{s))) RliO) = - [ i^^^i^R,(z)dzdz. 

Jc ^ 

We then proceed as in the Xi{s) case. 

To conclude the proof, we need to verify (I2.14p and (12.181) . However, we explicitly have 

[Xi(s), H{s)] = [P{s), gaiHis))]. 



15 



Thus, 



Y,{s) = [P{s), {P{s) - g^{H{s)))]. 



(2.26) 



Since P{s) + (1 - P{s)) = J, estimate (EUD follows from (^M) and (fT3|) . 
On the other hand, 



[Ms), H{s)] 



'l-ga{H{s))) R'M H{s) P{s) H{s) 



:i-gaiH{s))) RsiO) His) P(s) 



h. c. 



From Assumption [T] with P{s) smoothly nearly spectral and E{s) = 0, we have 

(His) P(s)) < 6/2. So, Assumption [31 and the bound f l2.17p imply that the norm of 

as 

le first contribution is O while the second contribution is equal to 

- {l-g,{H{s))) Rs{0) H{s) Pis) + 0{6/a) 



= - (l-gaiHis))) Pis) + 0(5 /a) . 
Putting everything together and using a < 1, we get 

[X,is),H{s)] = [Pis), gaiHis))] + 0{5/a') 
So, 1^2(5) defined by fl2.15p satisfies 

Y^is) = [Pis), iPis) - gaiHis))] + 0{S/a') 



(2.27) 



The rest of the argument is the same as for the l^i(s) case. 



2.2 Proof of Theorem D 

Existence of the propagator is a consequence of part (5) of Assumption [51 

We begin the proof of Theorem 11.21 by constructing a suitable family Pis). We note that 
one reasonable candidate for P(s) is a finite rank spectral projection Puis) for the Dirichlet 
restriction H^is) (extended by zero to the whole R'^): 

Puis) = - -i- / iHnis) - z)-' dz, (2.28) 

where T ■.= {ze([! : \z - Eis)\ = A/2 }. 
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Indeed, for x in the complement of dVL[s) and any vector 0, we have 

{H{s)Pn{s) <P){x) = E{s) {Pn{s) ^) (x) 

as desired. Unfortunately, the range of Pq{s) is not in the domain of H{s), since it is not 
twice differentiable at the boundary dQ{s) of the set Q{s). So, we must modify this family. 
The key estimate that will enable us to control the errors introduced by this modification 
process is encoded in Lemma [3.51 (c./. with the related result in |10j). 

To define the desired family of projections, let { ipi be an orthonormal basis for 
Range Pn(s), and define 

i'i = {I - Xs) ^i, 

where Xs is a smoothed characteristic function of the set dQ{s). By this we mean that Xs 
is twice differentiable as a function of x G JR!^, and if x G then Xs{x) = 1, and 

if dist{(9fi(s), x} > c/2, then Xs{x) = 0. We can assume Xs has been chosen so that 
||Xs||2,2 < where ||Xs||2,2 = / 1(1 - ^x)Xs{x)\^ dx. 

We define P{s) to be the orthogonal projection onto the span of { tpi }7=i- 

Lemma 2.3. Let 

Jn 

Q^'-Hs) := — (Pnis) - P{s)) . 

Then 

||g(")(s) II < Cn e-^/V^", (2.29) 
for n = 0, 1, 2. 

Proof Let R^{s) = (^{I - Xs) Pn{s) {I - Xs) - w ' and R^{s) = (Pf,(s) - ^^; /) \ 
for wE(C. Since (j(P(^(s)) = {0, 1}, Lemma 1331 shows that 

Pn{s) - (/ - Xs)Pds){I - Xs) = XsPnis) + Pnis) Xs - XsPn{s)xs = 0(6"^/^). 

Analytic perturbation theory then shows that for sufficiently small h and |ty — 1| = 1/2, we 
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have 

\\Rw{s) II < 3 (2.30) 



0(e-'^/^) (2.31) 



P{s) = - ^ RUs)dw and Pn{s) = " ^ RM dw , (2.32) 



where j := {w & (D : \w — 1\ = 1/2}. Combining these bounds, we get the estimate (I2.29p 
for the n = case. 

For the n = 1 case, we note that 

^R^{s) = - R^{s) Pn{s) Ru,{s); 
as 

^^Ru,{s) = - RUs) ^ { a - Xs) Pn{s) (J - X.) } RUs) . 
Using f l2.3ip . Lemma [3.51 and smoothness of Xs, we obtain that 

^ (^RUs) - RU^)) = 0(e-''/VA) , 

and the bound 02.291) for n = 1 follows from the contour representations (12.321) . The n = 2 
case is handled in the same way. H 



Proposition 2.4. Assume the hypotheses of Theorem M.^ and define P{s) as above. Then 
Assumption\J\ holds with 5 = C he~^/^. 

Proof It is clear from the definition of Pn{s) that (1 — Xs) {H{s) — E{s)) Pn{s) = 0. We 
therefore have 

{H{s)-Eis)) (l-Xs)Pnis) = - [H{s), Xs] Pn{s) = - [H {s) , F,] Pn{s) , 
where Fg = Xs ■ Xn{s) is a smooth function, supported in Q{s). Lemma [331 then shows that 

UHis)-E{s)){l-Xs)Pn{s){l-Xs)\\ < C h e-^'\ 
The result now follows from the identity 
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and the bound (^M)- 



Our next goal is to show that Assumption [2] is fulfilled. Our proof relies on the geometric 
resolvent identity (Lemma 13.31 in Appendix). 

Proposition 2.5. Assume the hypotheses of Theorem }!.^ and define P{s) as above. Then 
Assumption^ holds with S' = C h e''^^''' \\\ga\\\i/ ^, for a < A/2. 



Proof We first observe that if a < A, then ga{Hn{s) — E{s)) coincides with Pn{s). Since 



Pn{s) (1 - P{s)) P{s) = (P^(s)-P(s)) (1 - P{s)) P{s) = 0(6-"/'') 
by Lemma 12.31 we only need to verify that 

iga{H{s) - E{s)) - gaiHnis) - P(s))) (1 - P(s)) P(s) 

= igaiH{s) ~ Eis)) - gaiHnis) - Eis))) P(s) P(s) 



(2.33) 



is exponentially small. We note that by construction 



XBis) Pis) = 0, where B := {x e M"" : Xsix) = 1} H Qis) . 



So, we can multiply -P(s) in (I2.33P by Xn(s)\B{s) from the left for free. We now use Lemma 
13.31 with the choices = IR"^, Ai = fi(s)\P(s) and A = fi(s) with ^^.G supported on the 
set K C fi(s) with dist {K, Ai} > c'. By taking adjoints, this yields 



iHis) - z) ^ x^{s)\B(s) - iHnis) - z) ^ Xn{s)\B{s) 
= iHis) - z)-' [His), 0] iHnis) - z)'' • 



(2.34) 



Hence 



{iHis)-z) ^ Xnis)\Bis) - iHnis) ~z) ^ Pis) 



iHis)-z)-' [His),Q] iHnis) -z)-' Xnis)\Bis) Pis) 



< II iHis) - z)-' [His), Q] II XBis) iHnis) - z)-^ P(s) 



(2.35) 
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where in the last step we have used [H{s), O] = [H{s), 6] Xb{s) and Xn{s)\B{s) P{s) = P{s). 
We bound the first norm hj C h |Im2;|~^ with C that depends only on \z\, using Lemma [3.21 
To estimate the second norm, we bound 

XBis) iHais)-zr' Pis) 



< 



Xb{s) (Hais) - z) ^ Pn(s) 



By Lemma [2. 3[ the second contribution is bounded by C e ''/''/|A| for \z — E(s)\ < A/2. 
To estimate the first term, we compute 



^JiHnis) - z)-' Pnis)} - {(HUs) - z)-'} Pni-s) 

d 
ds 



-1} Pn{s) + {E{s) - z)-^ Pn{s) 



+ {E{s) - z)-^ (Hni-s) - z)-' Hni-s) Pni-s) . 
Note that Hn{s) = V{s) by Assumption [5l with suppy(s) C /(s). We can estimate 

Xb{s) (Hais) - zy'^ Pq{s) 



< 



\E{s)-z\^ 
1 



Xb{s) Pni'S) II + 



1 



\E{s)~z\ 



+ 



Xb{s) {Hn{s) -z) ^ XI 



\Eis)-z 

Assumption [5l and the first resolvent formula show that 



Xb{s) Pn{s 
V{s) Pn{s 



(2.36) 



V{s) (Hnis) - zy^ < C/A when \z ~ E{s)\ = A/2. 



< C. All other 



So, using the contour representation fl2.28p . we can bound V{s) Pn{s) 
terms in fl2.36p can be bounded using Lemma [331 (for \z — E{s)\ < A/2). We thus see that 



Xb{s) (Hnis) - z) ^ Pn{s 



< C 



A |Im^|2 

Putting everything together in (I2.35p . we obtain 

{{H{s)-zy^ Xn{s)\B{s) - {Hn{s)-zy^ Xm\B{s)} P{s) 

- |Im2|3 A ■ 

The lemma now follows from the Helffer-Sjostrand representation (12.201) . 



(2.37) 
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3 Appendix 



Here we collect a number of the technical statements used throughout the text. Many of these 
are well-known results that we have generalized to include magnetic fields and/or restricted 
domains. 

Let H ■.= 'Pa-Va + V on IR^, where Va = -ihV-A{q). Let Hq denote its Dirichlet 
restriction to the set Q. We assume H is self-adjoint and that Hq is bounded from below, 
so that it admits the Friedrichs extension. Let G and G be a pair of smoothed characteristic 
functions supported inside Q, and taking values in [0, 1], such that GG = G (which means 
that G is "fatter" than G). Throughout this Appendix, we assume that A and V satisfy 
the following hypotheses: 

Assumption 6. 

Al The components of the vector potential A and their first derivatives diA for i = 
1, d are bounded on IR'^. 

A2 For X in the support of VG, < C. 



We frequently estimate the norm of an operator A by looking at A* A. The following bound 
will be used throughout the Appendix: 

Lemma 3.1. Let Hq := 'Pa ■ 'Pa + 1 o^nd let Hq be either Hq or the Friedrichs 

-U- 

extension of its restriction to Vl. Similarly, let 'Pf denote either i-th component of "P a or 
its restriction to VL (which is not self-adjoint). If G is a smooth bounded function with 
support inside Vt, then for h < 1, we have 

^1/2 



Vt G (H* 



< \\G\\^ + Ch\\G\L. (3.1) 



Remark Note that in the lemma, T'j G = T'i G. 



Proof The non-trivial part of (13.11) is the bound for the Q restriction. (See e.g., [15] for 
the IR'^ case). The first step is to show that for any G^ function ( supported inside f2. 



2\-l 



< 



4 IICIIi.o. 
\t\ + 1 



(3.2) 
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for any value of the parameter t & M. To prove this, we write 

+ [H^ + eY' [C, Va ■ Va] C (^^o'' + t'Y' (3.3) 
The first term here is bounded by ||C||^ (1 + To bound the second term, note that 

[C^Va-Va]C = [CVa-Va] + [[C,Va-Va],C]- 

The second contribution, [ [(, 7^ a ■'Pa], C], is equal to 2 h"^ (VC)^- We use this to see that 
by making an error whose norm is bounded by 2 /i^ II Cll i oo /(-'- + ^^)^' second term in 
(]3.3p . which equals 



can be replaced by 



However, 



has norm bounded by (1 + 1^)^^ Putting the various pieces together yields (13.21) by some 
simple estimates and ^ < 1. 

Now we derive the bound (13. ip . First, we have 

2 



Vf G (H^) 



< 



< 



{H^)-'^' GVa-VaG {H^) 



We use the representation 



to see that 



[E^X'' G «) 



Q\-l/2 



G 
G 



[E^f' G,{H^) 



Jl\-l/2' 



(3.4) 



I WY^" I W + G] « + ey' dt. 
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Since 

[H^, G] = -2 i hV^ - (VG) + (A G) 
we can use the bound (13. 2 p with ( = djG to estimate 



K,G] (H^ + f)-' 



< 



Gh II Gil 



2,oo 



On the other hand, 



< 



\t\ + 1 

2 



|t| + 1 



so the right hand side of (13.41) is bounded in norm by || G 

Our next step is to estabhsh the following uniform bounds: 

Lemma 3.2. For the setup above, we have 
II [H, 6] {H-zy^ \\ < Gh(^l + 



oo + G^ II G II2 . 



Im2 



(3.5) 



with G depending only on 9. Furthermore, if E E M satisfies dist (cr(ifn), E) > A > 0, 
then 

1 + 1^1 



[Hn,e] {Hn-z)-'\\ < G h (^1 + ) 



(3.6) 



for any z & G , such that 1 2; — £"1 < A/2. Here G depends again only on O. 



Remark Note that if if) is smooth with support in fi, then [if, = [iJ^, G] %p. 

Proof Let denote either H or iff^, and similarly for the other operators that appear. 
We observe that 

Let Rf denote the resolvent in (13.51) (or in 13.61 accordingly) . then 



2 h 



\\[H#, 6] i?f II < fl'||A0|U ||i?f II 
To bound the second term here, we write 

V* ■ {ve) Rf = v*Q (H*y' H* ■ (ve) R* 



V* ■ (V 6) R* 



+ V*e (H*y' . [H*,VQ]Rf 



(3.7) 



(3.^ 



(3.9) 
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where Ho is defined as in Lemma 13.11 and G is defined as in the beginning of this Appendix. 
Note now that 



H* Rf = I + {{l + z)I - V#) Rf. 



(3.10) 



Using this, Assumption [HI and Lemma 13. H we see that the first term on the right hand side 
of fl3.9p is bounded by 



vi e (h* 



vi e (h* 



live 11,,^ + 

< C + C{l + \z\) \\R*\\ , 
where we have absorbed the Sobolev norms of G into C. 



{V&)-{{l + z)I-V#)\\ \\R*\ 



(3.11) 



To bound the second term in ( 13. 9p . we write 

= - (A 9,0) - 2 z n e e ■ (V diQ) . 

Using this, we see that 



V* e ( H* 



< c h' + c h 



vt e (h* 



2 h 



-1/2 



ivaei 



H*] e e 



< c h, 

where we have repeatedly used Lemma 13.11 with G = Q and the identity 

e (vty e = evte. 

We can consequently bound the second term in fl3.9p by 



V*e [H*y' ■ [V*-V*, V0] Rf 



< C h\\R* 



Putting everything together into (13. 8p . the bounds (13. 5p and (13. 6p follow from 

II II < l/IIm^l, and || i?^ || < 2/A whenever |E - ^1 < A/2. ■ 

Many of the technical results throughout the paper rely on the following general result 
(Lemma 4.2 in ^) that is known as the geometric resolvent identity. 
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Lemma 3.3. [The Geometric Resolvent Identity] Let H be a Schrddinger operator on 
IR!^. Consider four open sets Ai, A2, A, and that satisfy Ai C A, A2 C A, Ac fi, 
and dist{Ai U A2, A'^} > 0. Let Q be a smooth function which is identically 1 on a 
neighborhood of Ai U A2 and identically on a neighborhood of A'^. Given any restric- 
tions Hq and H\ of H to Q and A, respectively, we have 

XA, {Hn-zy' = XA, {Ha~z)-' e (3.12) 

+ XA, (H^-z)-' [H, 6] {Hn-z)-' 

for any z for which both resolvents exist. Also, 

XAi {B.^-zY^ XA2 = XAi {B.K-zy^ XA2 (3.13) 

+ XA, {H^-z)-' [H,Q] (Hn-z)-' xa, , 

under the same conditions. 

Proof Since the function B has support strictly contained within A, multiphcation 
by 6 satisfies Ha Q = Hq Q on V{Hn). Thus, on V{Hn), 

[H, e] = {HA-z)e - e {Hn - z) . 

We multiply this on the left by {Ha — z)~^ and on the right by {Hq^ — z)~^ to see that 

e {Hn - z)-' = {Ha - z)-^ e + {Ha- z)-' [H, 6] {Hn - z)-\ 

Multiplying both sides of the above equation by xAi from the left and using xAi ■ © = XAi 
gives (13.121) . Multiplying both sides of (13.121) by XA2 on the right gives (I3.13p . H 

Armed with this tool, we can prove 

Theorem 3.4. [The Gombes-Thomas estimate] Let H = 'Pa ■ 'Pa + V be as above. 
Suppose that on the domain J, the potential V is greater than E + b, for some h-independent 
b > 0. Then there exists an h-independent t] > 0, such that the resolvent of the Dirichlet 
restriction Hj of H to J satisfies 

II xj. {Hj - z)-' XJ, II < K (h-^ + e-^^l^ (3.14) 
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for z G (T, such that \E — z\ < b/2, and any Jij C J that satisfy 

dist {Ji, Jj} > c/2. 
Here K is a constant that depends only on c. 

Proof of Theorem 13.41 Consider the operator H := H^Va ■ a + V{x) acting on M'^, 
where 

' V{x) ii X e J 
E + b otherwise. 



V{x) 



Then the (improved) Combes- Thomas estimate [8] is apphcable for H and imphes 



XJ. [H-z 



(3.15) 



for \E — z\ < b/2 and dist {Ji{s), Jj{s)} > c/8, and some generic constant K. 

If we now use the geometric resolvent identity (13.131) . with XAi = XJi^ ^ = ^'^^ and 
A = J, where G satisfies 

dist {supp (Ve) , xjJ > c/8 , for i = 1, 2 , 



we get 

IIxji {Hj-zy^ x-h 



< 



X.h [H-z 



XJz 



(3.16) 



+ II xj, {Hj-z)-^ [H, 0] 
where J3 can be chosen in such a way that 

[H, 6] XJ3 = [H, 6] , and distjJg, J2} > c/8. 



XJs [H-z 



-1 



XJ2 



Using (13.61) and (13.151) . we bound the right hand side of (13.161) by 

I + \E\' 



K K 



.-1 



since the gap A in (13. 6p in the situation at hand is b. I 

The Combes-Thomas estimate leads to the following result, where we suppress the s- 
dependence whenever possible for the sake of brevity. 
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Lemma 3.5. [Exponential bounds for the spectral projection] Let Hq satisfy Assumptionl5[ 
Suppose F and its first and second derivatives are smooth and bounded, and that F vanishes 
on B := {xelR^ : dist {x, dVL] > c/8}. Then 



\\FPn\\ < Ce-''/^ 
\[Hn,F]Pnis)\\ < C h e-^/'\ 



and 



II F (Hn-z)-^ XiW < Ce-'^/^/A, for \z - E\ < A/2. 
Here I is a set that is slightly larger than I and satisfies I G I and 

dist {/,/"} > c, dist {S, 7^} > 3c/4 
for some h-independent c' > 0. The constant C depends only on c, c' and E. 



(3.17) 
(3.18) 
(3.19) 

(3.20) 



Proof of Lemma 13.51 Assertions fl3.17p - fl3.18p are consequences of (13.201) . so we prove 
( ]3.20p first. To do so, we make use of Lemma 13.31 We let A = f2f] J and Ao C A be a 
neighborhood of B. We choose G so that d^Q is supported in a set S* C J\Ii such that 
dist {5, B} > c/2. Lemma IX^ then yields 



F (Hn-z)-' 



Xi 



F{Ha-z)-' [H, 0] (Hn-z)-' 



XT- 



Since xs [H, Q] = [H, 6], we obtain 



F (Hn - z) 



XT 



{F (Ha-z)-' xs) {[H,e] (Hn-z)-' xt) 



(3.21) 



(3.22) 



Note that for 1^; — E'l < A/ 2, the first term in the parentheses on the right hand side enjoys 
the property (IXTD . Putting together flX^ . fIXTD . and ([XSD, we get ^?Mf . 
Note now that Xi H{s) = H{s), 



-1 



(Hnis) - zy' His) (Hnis) - z 



-1 



and 



H(s) (Hn(s) - z 



-1 



< C/A for \z-E{s) \ = A/2 
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by Assumption [51 It follows from the integral representation fl2.28p and the bound fl3.20p 
that the first assertion in flS.lSp holds. The second assertion in fl3.18p is established in the 
same way by taking the second derivative of the resolvent (with respect to the s variable). 

We now prove the bound fl3.17p . The integral representation fl2.28p and the bound fl3.20p 
show that 

II F Pn Xi II < C e-"/^ 
The desired bound fl3.17p will be obtained if we can show that 

II F Pn Xi II < C e-"^!^ implies \F P^\ < C e""/^ . 

To prove this assertion, it suffices to show that for every t\} G Range with = 1, we 
have 

IIX/^II > K > 0. (3.23) 
Indeed, it follows from f l3.23p that 
PnxjPn > K^Pn, 

so 

F Pn x\ PnF > {F Pn F) , 

and the result follows. 

Inequality ( I3.23P is a consequence of the bound 

||F,^|| > K, (3.24) 
for some smoothed characteristic function Fj of the set / that satisfies 

Fj Xi = Fj , and Fi{x) = 1 for x G / . 
To show that (13:241) holds, let 

ipi := Fj ip, := (/ - Fi) ip, and ip = ipi + ■ 

Then 

(^2, Hn^2) > (E + b) ||^2ir 
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by the definition of tlie region J. On tlie otlier liand, 



E 



< E 



(E + b) 



= E + 6||^2f + {[Hn,Fi]i;,i;,). 
We now use tlie integral representation f l2.28p together with the bound (13 .Gp to see that 



[H, Fj]iP\\ < C h II Fi 



2,00 



where C is independent of the choice of ip ^ Range Pq. So, with the appropriate choice of 
Fi, we see that || ^2 T < C h, and f l3.24p follows. 

Finally, we show that f l3.19p follows from f l3.17p . The proof closely follows the proof of 
Lemma 13.11 We use (13. 7p and (I3.17P to bound 



[Hn,F]Pn\\ < 2h 



VaFPu 



+ c e-^/^ 



where F := |V-F|. Now, we have 

PuFVa-VaFPu = Pn F^ Va -VAPn + Pn [[F, Va ■ Va], F] Pn . (3.25) 
The first term is bounded by 



Po F 



F Pn 



using the bound (I3.17p . the identity (I3.10p with the estimate in (13. lip , as well as the contour 
representation (ESSD- Since {F,'Pa ■ Va), F = 2h^ (Ve)^ the second term in i^TI^ 
is bounded by C e~^/^ using (I3.17p . H 
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